(N 



C/2 



m 
in 



X 



Metal-Insulator Transition of the Quasi-One Dimensional Luttinger Liquid 
Due to the Long-Range Character of the Coulomb Interaction 



V. S. Babichenko 

RSC " Kurchatov Institute", Moscow 123182, Russia, 
e-mail: babichen@kurm.polyn.kiae.su 



An instability of the quasi-lD Luttinger liquid associated with the metal - insulator transition is 
^ considered. It is shown that the homogeneous metal ground state of this liquid is unstable and the 

charge density wave arises in the system. The wave vector of this charge density wave has nonzero 
component both along the direction of the chains and in the perpendicular direction. The ground 
state has a dielectric gap at the Fermi surface, the value of this gap being calculated. 



1. The instability of a quasi-lD electron liquid with respect to the metal-insulator transition with the formation of 
1^ ■ the charge density wave (CDW) is analyzed in the present work. The materials, which can be considered as candidates 
5^ for the theory proposed below are quasi-lD organic conductors, for example, (TMTSF)2X in strong magnetic fields 
[1, 2]. 

The motion of particles is supposed to take place along the ID semiconductor chains which form the crystal lattice 
with the lattice constant A and represent a quasi-lD doped semiconductor. The thickness of the chains a is supposed 
to be much smaller than the distance between chains A (a << A ). The interaction between particles is the usual 3D 



' Coulomb interaction Vp^p — e*^/ | r — | with the static dielectric constant eo . So, the effective charge of particles 



e^/ eo where e is the bare electron charge. The analysis of the role of the long-range Coulomb interaction 
Q character is the main goal of the work. 
O For the first time the different types of instabilities of a quasi-lD Fermi liquid with the short-range bare interaction 

in the system of the ordered metallic chains in the parquet approximation have been investigated in the work [3] . The 
quasi- ID Luttinger liquid with the long-range Coulomb interaction has been investigated in the work [4]. In that work 
J> ■ only the forward scattering of particles has been taken into account. The influence of the long-range character of the 
Coulomb interaction in carbon nanotube which can be considered as consisting of two chains has been analyzed in [5] . 
In this case the Coulomb interaction between two chains does not change the Luttinger character of the ID electron 
liquid. The metal-insulator transition of the electron liquid in a strong magnetic field in the parquet approximation 
in the case when the Coulomb interaction can be considered as the short-range interaction has been considered in the 
^\ work [6]. Note that the long-range character of the Coulomb interaction can change the situation from the parquet, 
0^ when several channels have to be taken into account [10], to the situation with only one separated channel [9]. The 
long-range tail of the Coulomb interaction in the electron liquid in a strong magnetic field has been considered in 
[7]. In the present work the influence of the long-range character of the Coulomb interaction on the properties 
I of the backward scattering in the quasi-lD electron liquid is analyzed. In addition, the influence of the short-range 
. correlations is taken into account too. If the Coulomb interaction is supposed to have only the short range character, 
(~| ■ i.e. the electrons interact only at the same chain, the ground state has the properties of the Luttinger liquid. The 
O ' influence of the long-range part of the Coulomb interaction, i.e. when the interaction between electrons situated at 
the different chains is essential, results in the instability of the metallic homogeneous ground state of the Luttinger 
liquid. In addition, the ground state becomes dielectric and inhomogeneous. The charge-density wave (CDW) arises 
in the system. This CDW has the wavevector components both along the direction of chains and in the perpendicular 
direction. The wavevector component along the direction of chains equals kz = Q = 2pp/h where pp is the Fermi 
momentum pp = hnnX^ and n is the density of the Fermi liquid. The wave vector component perpendicular to the 
direction of chains equals | 1^ j_ \= \/2tt/\. In this case the amplitude of the CDW does not depend on the chain 
coordinate, but the phase changes its value from zero to tt for the neighboring chains. The ground state of the system 
is dielectric and has a dielectric gap at the Fermi surface. Below we use the system of units with h — I. 

The essential point of our analysis is the consideration of the interference between the channel of the backward 
scattering and the channel, which results in the infrared divergences, namely, the channel of the forward scattering. 

The transitions between different chains are neglected, but they do not play an essential role if they are small in 
comparison with the dielectric gap value arising at the Fermi surface. We suppose that each chain can be considered 
as a continuous medium and the system of chains forms the square lattice with the lattice constant A. The parameters 
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of the system is assumed to obey the inequahtics ^ ^ Pf ^ j where ub = l/m*e*'^ is the effective Bohr radius and 
m* is the electron effective mass. The first part of this inequality means that the non-dimensional constant a = e*'^ /vf 
, which defines the value of the effective interaction, is small a << 1 where vp is the Fermi velocity. Because of 
the condition a << 1 the random phase approximation (RPA) is correct for the calculation of the thermodynamic 
properties, for example, free energy [11]. The second part of the inequality means that the density of the electron 
liquid is sufficiently small so that the Coulomb correlations are large and the long-range character of the Coulomb 
interaction is pronounced. Below, for the simplicity of notations, the effective charge e* will be denoted by e. 

The electron liquid is supposed to be spinless. This situation can be realized by the switching-on the strong magnetic 
field directed along the chains. In this case the electron spins are frozen, but the motion of the electrons along the 
chains are not perturbed by the magnetic field. Moreover, the jumps between chains are suppressed by the magnetic 
field. 

2. The action of a quasi-lD Fermi liquid with the Coulomb interaction has the form 

S = Sq + Sint (1) 



Sint — — 2 2^ 



I dtdzdz' {i:{t,r)i;{t,r)) V?^?, {ip{t,r')'ip{t,r')) 



where z -is the direction of chains, Pz = —idz , r ^ ( z. It), n are the discrete coordinates of a chain in the plane 

perpendicular to the direction of chains. The field '4>{t,f) is the Fermi field (Grassmann variables). For simplicity, 
we suppose that the lattice is the square lattice with the lattice constant A. The motion of electrons from one chain 

to another is neglected. 

The correlation properties of the Fermi liquid are connected with the correlations of quasi-particles near the Fermi 
surface. The fields of these quasi-particles are denoted by V'^"'"^ for the part of the Fermi surface with = +Pf and 

for Pz = ~pf ■ 

The Coulomb interaction can be decoupled by the introduction of the virtual plasmon fields of and $ [7] . The 
plasmon field cf) has the small z-component of the momentum transfer kz and corresponds to the forward scattering of 

quasi-particles and the plasmon field <I> has the z-component of the momentum transfer kz close to the value Q = 2pF 
and corresponds to the backward scattering. The action of quasi-particles can be represented in the form 

Sl[^, <t>, = Sf[^, <P, + [^] + sill m > (2) 

where 



(3) 





SfLard {<l>]=-\jdtj dkzd^k^cj>{t, -kz, -k^)V-' (kz, kj_) Cj>{t, k„k^) (4) 




/3 

silk = " / '^^^ / dkzd^k^^^{t, kz,k^)U-^ {k^_) kz,kA_) (5) 



In these expressions the momentum k ± belongs to the elementary cell of the reciprocal lattice formed by the system 
of chains. A; _l is perpendicular to the chain direction z and /? is the inverse temperature. The temperature is supposed 



2 



equal to zero, but the Matsubara technique is used. The functions G^^-* are the Green functions of quasi-particles 

near the different parts of the Fermi surface with pz = ±Pf in the external plasmon field (j) , Pz = —idz and ^ are 
the discrete coordinates of chains in the plane perpendicular to the chain direction. 

(Of^y' =-drTVFPz-i(t>{t,Z,^), (6) 

The function V ^fc^, k±j is the Fourier transformations of the Coulomb potential with the small momentum transfer 

kz , and U (ji±^ is the Coulomb potential with the component of the momentum transfer kz close to the value of 
kz=Q = '^Pf 
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Here 6 are the vectors of the reciprocal lattice so that b = (27rn/A; 27rm/A) where n, m are the integer numbers, and 
fc_L belongs to the cell of the reciprocal lattice (— tt/A < k^- y < tt/A) . It is convenient to extract the term with 6=0 
in the sums (7) . As a result, the sums (7) can be represented in the form of sums of two items. The first one is the 

summand with "6* = and the second is a sum over ^ in (7). Thus U {^A^j and V (kz^, ^A^j are represented in 
the form of the pseudo-potentials having long-range and short-range parts 

v{kz,kA)= -^^^ + Vcore {kA_) (8) 
U (kl) = + C/core (9) 

The functions Vcore {k±j and Ucore (^^) can be considered as independent of fcj. if | k±_ \« A. In the case of the 

small thickness of chains a « \, which is supposed, the dependence of Vcore and Ucore on the momentum fcj^ is weak 
for I fcj^ |~ A too and thus can be neglected 

Vcore (fc±) « Ucore (fc±) ~ Kore (0) = Ucore (0) = ^ (10) 

The value A in (10) has the form 

P = (2/7r)AMn-i , 

where A = tt/A, Aqo ~ l/o, thus Aqo >> A and A << A. The parameter Aqo is the cutoff momentum in the sum 
over 6 in (7). The truncation of the sum (7) is the necessary operation in the connection with the divergency of 
this sum at the large momentums b . The Fermi momentum is supposed to be siifficicntly small, so the inequality 
Q « A << A obeys. The functions V {kz^k^^ and U {^a^ (7) as well as V^'^ [k-z^kA^ and J7~^ (^^) ^^"^ ^'^'^ 
periodic functions of k± with the periods 1} x,y = (27r/A; 0) ; (0; 27r/A) and k±_ in the expressions (8), (9) belongs to 
the elementary cell of the reciprocal lattice. The function U~^ (^-'-) represented in the form 



C/-1 (k,) = U-,lc (k^) 



t^+Q'^+Aire^Ucole [k_ 



(11) 



The coordinate representation oiU ^ (^-'-) ^^'^ written as 
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where C j is 



The region of integration in (13) is the elementary cell of the reciprocal lattice. The value Aire^l/'gl^ (^-'-) ^^^isfies 
the estimations Attc'^U^J^^ (^-'-) ~ ^Tre^C/^J^^g (0) = << . Thus, the momentums k±, which give the main 
contribution to the integral (13) for the function C i are of the order of A << A. For this reason, the function 
^core (^-l) equality (13) can be taken equal to f/j^J-e (0) . Thus, we obtain 

C (i?) = ^Ko (1 ^ I Aq) for I i? I a » 1 (14) 



where Kq (x) is the cylinder function of imaginary argument and Ag = \J A'^ + Q"^ . The value of | ^1 — ^2 | in the 
equality (12) is supposed to obey the inequality | ^1 — ^2 |>> 1/A. If | ^1 — ^2 |>> l/Ag ~ 1/A , the function 
Kq ^1 — ^2 I Agj can be replaced by 

Ko{\ t,-%\lQ)^j^^6 
In this case the equality (12) can be written as 

(''"")7?.,7?.-5^(o^-^bX^'-S0 

The expression (15) is correct in the case of the slow change of the fields $ at the scale of the order of the distance 
between neighboring chains. When the characteristic distances | ^1 — ^2 | satisfy the inequality 1/A <<| — 
^2 \« 1/A , the non-local character of the function C ^| ^1 — ^2 |^ in (12) is essential. Using the equality (12), 
the expression for (^~^)^ ^ can be represented in the form 

The general property of {U^^)-^ is the negative sign of (t^^^)^ ^ ^1 7^ ^2 under the condition QX « 1. 
3. Green functions g'"^^ can be calculated exactly [8] 

Gf (x, y, H) = exp (^9^^^ (x) - (y)) G^^^ (x - y,lt) , (17) 
where x = {t, z) and the Fourier transformations of ^'^^ (x) are connected with the plasmon field (f) by the equation 

In the Eq. (17) the functions G^^^ (^x — y, are the Green functions of the non-interacting quasi-particles located 
at the chain 
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iuJTVF (Pz TPf) 



Our goal is the calculation of the effective action for the plasmon field To obtain the effective action we consider 
the generation functional Z 



Z = j Di;D^D<j)D<^+D<S> exp {Sl + Ssources} , 
where Sgources depends on the sources J, J and has the form 



(19) 



Ssources (fx $+ (x, ^) J (x, itj + J (^X, ^) $ (x, ^) 

The integration over the Fermi fields in the expression for the generation potential Eq. (19) with the use of the Eq. 
(2) gives 



Z = j Di)Di)D(i)D<^+ exp {Spi + Ssources} : 



where 



spi = snet [0, + [-A] + siz, m + 



(0) 



sources 



(20) 



(21) 



Soet [4>, *f ] = SpLn 



(22) 



Note that for $ = the part of the action Suet \4>i ^] can be calculated exactly [8]. As a result, the action Spi can be 
transformed to the form 



where 



Spi = Soet [A] + Sf 

orwara 



sources 1 



Snet [A] = SpLr 



(gS*')" 



-zA 



(23) 



(24) 



Sfo 



I 



(kodkzd k± 



6{-u),-kz,-k±)DQ ^ (u),kz,k±^ 0{u),kz,k±) 



(25) 



S^-^e,A] = - [d'x Y: e-i'"'^)A+{x,^)u-'(%^')A{x,^')ei''''^') 

In these expressions the new fields A , A+ are introduced instead of the fields $ , They obey the equalities 
a(x,;S) =$(a;,i?)e"^(''''^) ; A+ (x,^^ = ^+ (x,Tl^ ei'''^ ■ 



(26) 



The field 6 is connected with the fields ^^^^ by the following way 



(k, u, i?) = 6'(+) (k, u, i?) - 0^-'> (fc, w, i?) = -2i 



vpk 
ijp' + {vpk)'^ 



(27) 



(28) 
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The inverse propagator of the ^- field Dq^ (cj,kz,k±^ has the form 
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£>o U,fcz,fc±) = ^ , , ,2 V-Ukz,kA ~ILi^\uj,kz) (29) 

^ ^ 8{VFkz) L \ / J 

where IIq'^^ ^w, fc^, is the zero polarization operator with the small momentum transfer kz {kz « Q) 

where is the squared screening momentum ( = 4e^/A^WF = 4aA^/7r^). 

Note that the field (t, z, has only the real values. This statement is a simple consequence of the expression 

(28). The expressions (24) and (26) can be obtained, for example, from the consideration of the diagrams which are 
the terms of a series expansion of (22) in powers of the plasmon field 

The simple expansion of S oet [A] in powers of the field A up to second power alone gives 

Soet [A] « A+HqA (31) 

where Hq = - {K^/2tt^vf) In {ep/ I A |) = - {t^ /Sne^) In {ep/ \ A |). 

4. First of all, we consider the approximation for which the dependence of Vcore ^^'^ Ucore {^'^ ((8)j (9)) 

on k±_ is neglected. This dependence will be considered below in the part 5 of this paper. In this approximation the 
values Vcore and Ucore Can be considered as equal values (due to the condition Q\ << 1 ). This allows us in the case of 
the spinless liquid to neglect Vcore ^md Ucore a-nd to take into account the long-range part of the Coulomb interaction 
only (first terms in the sums (8), (9)). This approximation is absolutely correct in the case of small momentums 
fc_LA << 1. In this part of the paper we suppose that this approximation can be used for /;;j_A ~ 1 too. Thus the 
action Eq. (23) can be represented in the form 

S^'r^ [6, A] = Snet [A] + Sf^:l, [6] + 5^17^ [A] + S^r^ [9, A] (32) 

where Suet [A] has the form Eq. (24). The terms S^forwardl^] ^'^'^ ^bt'ck ^ [^] have the forms (4) and (5), corre- 
spondingly, in which V~^ (^z,k±^ and U~^ (^-'-) ^'^o^^'i substituted by V~^ (jtz,k±j = (kz'^^'j^ /47re^ and 
U^^ (k-ij = + /47re^. Besides, the field A has to be substituted instead of $. 

In this case the part of the plasmon action s[^^^'^ [9, A] , which describes the interaction between the fields 9 and 
A, can be represented as a sum of two terms 

where 

i2 



ST' [^^ A] = dtdz [a+a {y^ey 



(34) 



5^* [e, A] = i A_ ^ j dtdzd'Tl (J (^±0)) (35) 

t = -z (a+ (^^a)-A(^^A+)) 

Note, that the presence of the term S*™* [9, A] leads to the divergency of the functional integral for the large values 
of '^_l6' and, for this reason, this means the necessity of the consideration of the large momentum k ± contribution 
or the lattice formulation of the plasmon action (23). The lattice formulation will be given below in part 5. 
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Thus the effective action for the plasmon fields can be represented as 

4f "')[A, 9] = Soet [A] - j eD-^e 



j dtdz A+ (q^ - ^1 - {^i.ef^ A 



+ 



d^'i^ [(a^ (^.a) - A (^,A+)) (V.O) 



Here Soet [A] is determined by the equality (24). It is convenient to introduce the amplitude and the phase of the 
field A in the following way 



A = pe'^ 
A+ = pe-'^ 



(36) 



Thus we obtain 



s^io^ = SpLn [g-1 (A)] - eD^'e -^Y.j {q^p" + (y^pf - (^^ + ^x))' 

where £>□ ^ is determined by the equality (29) and Soet [A] = SpLn G^^ (A) (Eq. (24)). 
We introduce a new variable instead of Q 



(37) 



r] = e + ix 



(38) 



After that we have 



(39) 



where the operator acts on the fields which are situated at the right side of this operator and acts on the 
fields which are situated at the left side of this operator. 
The integration over 77 gives 



S,io^ [A] = SpLn [g-1 (A)] - ]^SpLn [D-^ (p)] - -^Y, j (qV' + 



(40) 



where 



The last term in (40) in the approximation of the small value of p^ can be represented in the form 



(41) 



After that we obtain the following expression for the effective plasmon action 



7 



S,io^ [A] = Spin [g-' (A)] - ^SpLn (p)] - -^Y, j ^tdz (q^p^ + (^^p)' 



(42) 



47re2 



There is nonzero saddle-point field for this action. The saddle-point equations in the case of the constant p as a 
function of t, z and ^ coordinates and for the field x have the form 



1 f dud^d 
Vf J (27r) 



+ 



1 r dujd^d'^k^ y-^ 



47re2 



47re2 



(43) 



= (44) 
The possible solution of these saddle-point equations can be represented as 

P = ipo 

where po is the constant real number field and the solution for the phase % is 

X = (45) 
Here is an arbitrary vector lying in the cell of the reciprocal lattice. The value of po has the form 

{Q'+tiy 



Po = ep exp 



2a 



(46) 



The statistical sum Z has the maximum value for 1^ ± = {-k/\,%/\). Thus the ground state of the system is 
inhomogeneous with a charge-density wave which has the wave vector along the direction of chains equal to 2pF and 
in the perpendicular direction equal to ll ±_ = (tt/A, tt/A). 

The fluctuations 5p of the field p can be taken into account by the following representation of the field p 

p = ipo + 5p 

Note that in the approximation of the slow variation of the fields p and x in the time t and z-coordinate the 
calculation of the first two terms in (42) gives 



SpLn G-i(A) 



-e/ 



dtdz 



Po 



SpLn[D-Hp)]=^^El 



dtdz 



Po 



5. In this part of the paper the dependence of Vcore [J^-^j '^"^ Ucore '-'^ ^J- ^''^ taken into accoimt. This 

is essentially for involving of the short-range fluctuations. If there is an instability of the system in the plasmon 
channel with the momentum transfer Q = 2pF , the plasmon effective action, obtained after the integration over the 
field 9 and dependent on the field A alone, tends to zero near the instability point. Thus, near the instability point 
the exponent in the expression (20) can be expanded into a series in action Soet [A] -|- 5*"* [0, A], and after that the 
integral over the field 6 can be calculated due to the Gaussian form of this integral. In this case the contribution of 
gznt pg^jj represented in the form 



< 5'"* [9, A] >e-- 



fd^x Y (2^' ^) i?') A (x, Tt') Ce [st) 



(47) 
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where 



--< e 



>e 



Here 5^ = ^, 



< e 

The integration over 6 in (48) yields 



>0= j DOexp (a;, 'M') - (x, ll^ + S forward [0]} 



>0= exp • 



1 f dujdkz(Pk± 



(l - cos (lcj_6'^^^ D{uj, kz, k±) 



The propagator D{lj, k^, kj_) obeys the Eq. (29). The calculation of the integral (49) for Co (^^^^ gives 



The function / has the form 



< e V /V / >0= exp 



f{5Tt) =0 



K 



(« I Sl^ l) 



for I S'l \» 1/Aoo 

for 51? = 



(48) 



(49) 



(50) 



(51) 



7 = 2/7r; ^ = In (£f/ | A |) ; Aoo = 1/a. 

Thus, the part of the effective action < S'*"* [9, A] >$ leads to the renormalization of ^ j in (23) and the 

renormalized quantity f/^^j ^ ^ has the form 



U. 



(52) 



Thus 



The calculation of the Fourier component of U^^^j, gives 



^e// (^) = 4^e2 



1-A 



A^ \ / A \ 



A, 



Q . 



A. 



where 



A/Aq 



A (kl, i)= j dxJo (x\kl\ /Aq) Ko (x) 



7a{+l 



(53) 



(54) 



(55) 



The momentum Aq is Ag = y A^ + Q2 Assuming A/Kq << 1, the above limit of the integral (55) can be replaced by 
the infinity. Note that i/^^j decreases with the increase of the value of | fci | for the small | Fl | (| fc^L | /Aq << l). 
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The definition of A (jt±, (55) gives A ^fci = 0, ^ = 0^ =1. Thus in the approximation << 1 and | k± \« 1/A 

we can put A ^Fl,^^ equal to unity A ^fci,^^ = 1 in (54). 

The effective action of the plasmon field A can be represented in the form 



Seff = SpLn G-i (A) + < 5™* [e, A] >8 



(56) 



< 



5-* [9, A] >e= -jd^^j kl) Kf) {^1) A (x, kl) 



Expanding SpLn 



G ^ (A) in the effective action Sgff (56) in field A, we obtain the quadratic action for the A field 

(57) 



and the inverse propagator of this field has the form 

r-^ (kl, ^) = -Hq + u-f) (kl, ? 

Here Ilg = (— k^/Sttc^) ^ where = iaA'^/ir^ and ^ is the logarithmic variable ^ = ln(eF/ | |) • 

Now we will consider the pole of the propagator F when fc _l = 0. Note that the situation with A;_l ^ is analogous 
to that which has been considered in the part 4. For fcj, = the equation determining the pole = In {cp / \ ojq \) 

of F (t^ = 0, e) 

takes the form 



-47re^nQ + 



1 




= 



(58) 



Using the expression for FEq and = (2/7r)A^ In (^) the Eq. (58) can be represented in the form 

2 / 



l70«ln(%l + 







(59) 



It can readily be seen that the right-hand side of the equation (58) is positive for ^ = 0. However, as ^ increases, the 
right-hand side of (58) becomes negative. Thus, the solution of this equation exists and for the supposition Q << A 
has the following value within the logarithmic accuracy 



(Jo |~ e^exp 



2Q2 



A'^'ja In ( 




ep exp 



2a 



(60) 



Note, that the Eq. (60) is correct when ^2§~ ~ Q^/k^ << 1- The solution (60) coincides with the solution (46) for 
^± = 0. Here loq corresponds to the pole of F = 0, The existence of this pole means the instability of the 
system. As a result of this instability there is a CDW in the ground state of the system. The z-component of the wave 
vector of this CDW equals 2pF . As it has been discussed in part 4 the wave vector of this CDW has the component 
perpendicular to the direction of the chains and equal to = (tt/A, tt/ A). The existence of the CDW results in the 
existence of the dielectric gap at the Fermi surface. The value of the gap is equal to A^ = po, where po is defined by 
(46). 
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